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ABSTRACT 

We combine published optical and near-infrared photometry to identify new low-mass candidate members in an area of about 0.64 deg^ 
in Corona Australis with the 5 -parameter method. Five new candidate members of the region are selected. They have estimated 
ages between 3 and 15 Myr and masses between 0.05 and 0.15 Mq. With Spitzer photometry we confirm that these objects are not 
surrounded by optically thick disks. However, one of them is found to display excess at 24 yum, thus suggesting it harbors a disk with 
an inner hole. With an estimated mass of 0.07 Mq according to the SED fitting, this is one of the lowest-mass objects reported to 
possess a transitional disk. 

Including these new members, the fraction of disks is about 50% among the total Corona Australis population selected by the same 
criteria, lower than the 70% fraction reported earlier for this region. Even so, we find a ratio of transitional to primordial disks (45%) 
very similar to the value derived by other authors. This ratio is higher than for solar-type stars (5-10%), suggesting that disk evolution 
is faster in the latter, and/or that the "transitional disk" stage is not such a short-lived step for very low-mass objects. However, this 
impression needs to be confirmed with better statistics. 

Key words. stars;low-mass, brown dwarfs - stars: formation - stars: pre-main sequence - stars: luminosity function, mass function 



1. Introduction 

In the last two decades, our understanding of the formation and 
early evolution of low-mass stars and substellar objects has been 
remarkably improved, but there are still many unknowns regard- 
ing the exact mechanisms at play in these processes (see e.g. 
Whitworth et al. l2007l for a review). The census of the low-mass 
population for an increasing number of star-forming regions is 
important to address the problem of the universality of the ini- 
tial mass function (IMF). The comparison of the disk fractions 
(that is, the number of sources with disks with respect to the total 
number of members) in clusters of similar and different ages also 
helps to constrain the timescales of disk evolution (e.g. Haisch 
et al. 1200 11 1. Moreover, the disk fraction and the spatial distribu- 
tion of stars with and without disks provide important clues to 
understand the star-forming history and the disk evolution of a 
particular region. 

Most of our current knowledge about these issues comes 
from the study of star-forming regions and young clusters at dis- 
tances within 100-300 pc from the Sun (see the reviews collected 
in Reipurth 2008). However, in the vast majority of these studies 
the identification of cluster members is based on the detection 
of features related to accretion and circum(sub)stellar disks (e.g. 
Ha emission, infrared excess), or indicative of strong activity 
(e.g. hard X-ray emission). Thus, these surveys are likely to be 
missing a fraction of the young low-mass population, especially 
objects that have dissipated all or most of their disks, that are 
in quiescence at the moment of the observations, and/ or whose 



activity level lies below the detection threshold of current X-ray 
surveys. 

In a recent paper Comeron, Spezzi & Lopez Marti (2009, 
hereafter CSL09 ) reported on a large-scale optical survey of the 
Lupus clouds, selected to include most of the areas surveyed in 
the Spitzer Legacy program "From molecular cores to planet- 
forming disks" (c2d). Combining their data with 2MASS near- 
infrared photometry, these authors developed an analysis proce- 
dure based on a dimensionless parameter S that allowed them 
to easily identify very low-mass members of the dark clouds as 
an excess in the expected S-parameter distribution of field cool 
dwarfs and giants (see Sect.[3]for details). In this way they were 
able to unveil a large population of objects belonging to Lupus 
I and III, which seems to be composed of very low-mass stars 
and brown dwarfs that have lost their inner disks on a timescale 
of a few Myr. The discovery of a substantial and even dominant 
population of thus far unnoticed members of one of our nearest 
star-forming regions stresses the important unknowns that still 
subsist in the observational characterization of young very low- 
mass objects. It also poses the important question of whether 
similar populations of yet unknown very low-mass members ex- 
ist in other regions. 

In this paper, we apply the S -parameter analysis to a set of 
optical and near-infrared observations of Corona Australis, an- 
other of the nearest regions with ongoing or recent intermediate- 
and low-mass star formation. Our aim is the identification of new 
very low-mass candidate members of this dark cloud. We focus 
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on the core in the direction of the intermediate-mass star R CrA, 
which contains a compact embedded stellar cluster known as the 
"R CrA cluster" or "the Coronet cluster" (Taylor & Storey 1984|. 
With an age of about 3 Myr, its distance has been estimated 
to be within 50-170 pc (see Neuhauser & Forbrich ,2008, for a 
review). The young stellar population in this region has been 
studied with different techniques, including Ho- surveys (e.g. 
Marraco & Ryndgren 1981), optical spectroscopy (e.g. Walter 
et al. 11997b . near-infrared mapping (e.g. Wilking et al. [1992l l. 
mid-infrared surveys (e.g. Olofsson et al. ll999b and X-ray obser- 
vations (e.g. Neuhauser et al. 120001 Forbrich & Preibisch 2007). 
Some brown dwarfs and brown dwarf candidates have been re- 
ported (Wilking et al. 1997, Fernandez & Comeron 2001 , Bouy 
et al. 120041 Lopez Marti et al. I2005I I. Very recently, this cloud 
has also been the target of Sp itzer observatio ns (SiciUa-Aguilar 
et al. 120081 Lopez Marti et al. in preparation| l. 

The structure of the paper is as follows: In Section |2] we 
present the data used in our analysis. Section [3] summarizes the 
basis of the 5 -parameter formalism, discusses the contamina- 
tion expected in our sample and describes the object selection. 
Section[4]is a discussion of the membership status and the prop- 
erties of our new candidate members. Finally, in Section [5] we 
draw our conclusions. 



2. Optical and near-infrared data 

The analysis presented here is based on the same set of opti- 
cal observations of the Corona Australis star-forming region re- 
ported in Lopez Marti et al. (i2005i . The survey was performed 
with the WFI mosaic camera at the MPG/ESO 2.2m telescope 
of La Silla Observatory. It covered an area of about 0.64 deg^, 
distributed in two fields of 34'x34', which are placed in the core 
of the Coronet cluster and on a region of lower cloud density to 
the South of it (Fig.[T]i. For the present work, only the broadband 
observations in the R and / filters were considered, which cover 
the dynamical ranges 11 < R <22 and 12 < / < 23. The survey 
is complete down to /? ^ 20 mag and 7^19 mag. For details on 
the observations, data reduction and calibration, see Lopez Marti 
et al. (i2005 ) and references therein. 

The optical data were cross-matched with JHKs photome- 
try from the 2MASS archive (Skrutskie at al. 2006"), which is 
complete down to 7 ^ 16.5. For the present analysis, we require 
that the objects are detected in all five RIJHKs bands. The to- 
tal number of sources fulfilling this criterion amounts to 6136. 
Our sample is complete down to a mass of about O.O15M0 at a 
distance of 130 pc and an age of 3 Myr, according to the Lyon 
models (Baraffe et aL |T998l Chabrier et aL |2000t . 



3. Analysis 

3.1. SED fitting 

The spectral energy distributions (SEDs) of the objects in 
our optical-NIR catalog were fitted to the model spectra of 
H auschild t et al. ( il999) and A llard e t al. (12000) as explained 
in lCSL09l (see also Spezzi et al. l2007l l. The observed magnitude 
at each band is expressed as 
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Fig. 1. Spatial location of the previously known candidate members 
(triangles) and the new candidates from this work (circles) within the 
Corona Australis dark cloud. Note that only those previous candidate 
members selected with the same criteria as our new candidates are plot- 
ted. The two WFI fields are also indicated. The objects are clustered 
around the intermediate-mass stars, which suggests that they belong to 
the Coronet cluster. 



ature; and S is a wavelength-independent scaling factor defined 

as 



5 = 5 log D{pc) - 5 log/?(/?o) - 5. 



(2) 



where M^^.,,, is a synthetic magnitude in the same passband, com- 
puted from the model spectra, for a \Rq star at the same temper- 



The ratio Ay/A^ is a constant given by the adopted extinction 
law (Cardelli et al. 1989, with Ry = 3.1). A value of log ^ = 4.0, 
thought of being representative for young low -mass objects, was 
assumed for all the sources. 

Assuming a certain T^ff for each star, an equation of the 
form ([T]i was set for each passband and solved by least-squares 
to derive Ay and S . The temperature that minimized the residual 
of the fit was then taken as the temperature of the object. The 
radius R was computed from the corresponding value of S for 
each source using Eq. Q and a distance value of 130 + 20 pc to 
the Corona Australis star-forming region (Marraco & Ryndgren 
'1981 ; Neuhauser & Forbrich 120081 1. The luminosity L was then 
computed with the relation L = 4no-R^T^jj., where cr is the 
Stefan-Boltzmann constant. Finally, the best-fitting Te/f and the 
derived luminosity were used to estimate stellar masses and 
ages from the evolutionary models and tracks by Baraffe et al. 
(11998) and Chabrier et al. (I2000I I. again assuming a distance of 
130 ±20 pc. 

We note that the use of a different distance would modify 
the estimated radii, luminosities and ages, but not the derived 
effective temperatures, extinctions and 5 -parameter values. For 
instance, a distance of 170 pc, as suggested by Knude & H0g 
(II998I I. would imply higher luminosities (by 0.23 in logarithmic 
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Spectral type 

Fig. 2. Tfff vs. spectral type for the previously known members and 
candidate members of the Corona Australis star-forming region listed 
m Table|3]T] Filled and open circles are used for spectral types derived 
from spectra (Sicilia-Aguilar et al. I2008t and from narrow-band pho- 
tometry (Lopez Marti et al. I2005t , respectively. The solid line is the 
temperature scale for PMS objects from Luhman l\999S and Luhman 
et al. J2003b . The dashed line is this scale shifted 200 K towards lower 

Teff- 



scale), larger radii (by less than 0.5 for R < 1.5Rq) and younger 
ages by about 2.5 Myr. However, because the ages derived using 
only photometry and isochrones, as in our case, intrinsically con- 
tain an uncertainty of a few Myr (e.g. Baraffe et al. 2002 2009 ), 
the distance uncertainty is not a major issue for our parameter 
estimates. 

The residuals for the best fit give in turn an estimate of the 
quality of the solution. We found that 94% of the objects have 
residuals below 0.2 mag, which is considered acceptable given 
the accuracy of the photometry. The robustness of the solution is 
supported by the distribution of values, which shows a clear 
minimum at the best-fitting value for each parameter, and a steep 
rise after the three best fits, especially for Tgjf. The differences 
in the residuals between the three best fits are < 0.01 in most 
cases. 

As discussed by ICSL091 poor fits will be obtained for unre- 
solved pairs and for variable stars (because the WFI and 2MASS 
observations are not simultaneous). In addition, the spectral en- 
ergy distributions of some objects may be affected by the sig- 
natures of accretion and warm circumstellar dust in the visible 
(veiling and strong emission lines) and near-infrared (infrared 
excess), leading to an erroneous estimate of the temperature. In 
both situations, our fits tend to overestimate both the temperature 
and the extinction (see !CSL09i for details). 

The mean extinction error was AAy =0.15 mag according to 
the outcome of the fitting procedure, suggesting that most ob- 
jects in our sample are located in the foreground rather than in 
the background of the dark cloud. Sicilia-Aguilar et al. (I2008I I 
provide independent extinction measurements for some cloud 
members, which are derived from their spectral types and near- 
infrared colors. For the sources in common with that work, our 
values of Av tend to be higher than those provided by Sicilia- 
Aguilar et al. However, the difference is not larger than 0.5 mag 
for most of the objects and can be attributed to variability and/or 
to the use of a different extinction law. 

The uncertainty in the temperature determination is esti- 
mated from the comparison of the three best-fitting results for 
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Fig. 3. (Tgff, S ) diagram for those sources in our survey with a derived 
effective temperature below 4500 K. Mass tracks and isochrones from 
the models of Baraffe et al. ( il998t are also plotted. 



each object, which generally span a range of 200 K in effective 
temperature. We thus take 200 K as the error in our temperature 
values. This agrees with the results of Spezzi et al. (|2Q07| and 
ICSL09I who estimated an error of 250 K on average from the 
comparison of their fitting results (obtained with the same pro- 
cedure) with the effective temperatures derived from spectra for 
a sample of well-known members of Chamaeleon II and Lupus, 
respectively. 

In Fig. |2] we compare our 7^// values with the published 
spectral types, derived from narrow-band photometry (Lopez 
Mart! et al. I2005I I or spectroscopy (Sicilia-Aguilar et al. 2008) 
for a sample of previously known members and candidate mem- 
bers of Corona Australis, selected in Sect. |3.3| below and listed in 
Table |3.1| We also plotted the temperature scale from Luhman 
(1999) and its extension to late M-type objects by Luhman et 
al. (2003). As seen in this figure, most of the objects are sys- 
tematically cooler according to our fitting than the predictions 
of the Luhman scale. Indeed, other authors have suggested that 
the true temperatures of PMS stars are lower than the predictions 
of this temperature scale (e.g. Mohanty & Basri 2003; Barrado 
y Navascues et al. 120041 ). the differences amounting to 150- 
200 K. Indeed, we find that shifting the Luhman scale 200 K 
towards lower temperatures (dashed line in Fig. |2]) provides a 
better agreement with our derived r^yy values. This difference is 
still within our estimated error range. 

3.2. Contamination 

In Fig. |3] we show the [Teff, S) diagram for the sources in 
our survey. We also plotted two isochrones from the models of 
Baraffe et al. ( 1998 ), corresponding to the expected location of 
objects of 1 and 20 Myr, and several mass tracks from the same 
models for masses between 1.4 and O.O2M0. This diagram sug- 
gests that no Corona Australis members are present in our sur- 
vey with temperatures between 4100 and 3400 K, as this range 
of Tgff is devoid of objects in the area between and above the 
model isochrones. Indeed, cloud members with T^jj- > 3400 K 
(spectral types earlier than about M3) are generally expected to 
be saturated in our optical survey (/ < 12 for a cluster age of 
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Table 1. SED fitting results for the previously known candidate members of Corona Australis selected with our analysis 



Name 


Ay (mag) 


Teff (K) 


5 




R(Ra) 


M(Mo) 


age (Myr) 


CrA 205 


0.04±0.14 


3000 


7.69+0.07 


0.0103+0.0006 


0.376+0.012 


0.072 


8.0 


CrA 452 


1.14±0.11 


3400 


6.07+0.05 


0.0755+0.004 


0.79+0.02 


0.4 


12.7 


CrA 453 


2.00±0.12 


3200 


7.39+0.06 


0.0177+0.001 


0.433+0.012 


0.11 


8.0 


CrA 465 


1.42±0.26 


2800 


7.89+0.13 


0.0065+0.0008 


0.34+0.02 


0.03 


1.1 


CrA 468 


0.12±0.15 


3100 


6.98+0.07 


0.0226+0.0015 


0.52+0.02 


0.13 


8.0 


CrA 4107 


0.73±0.15 


3000 


6.63+0.08 


0.028+0.002 


0.62+0.02 


0.072 


2.8 


CrA 4108 


0.16±0.16 


2900 


6.80+0.08 


0.0205+0.0015 


0.57+0.02 


0.055 


2.2 


CrA 4109 


0.45±0.13 


3000 


6.24+0.06 


0.039+0.002 


0.74+0.02 


0.075 


1.4 


CrA 41 10 


0.94±0.16 


3000 


7.09+0.08 


0.0170+0.0013 


0.50+0.02 


0.07 


4.0 


CrA 41 11 


0.22±0.22 


2900 


7.50+0.11 


0.0107+0.0011 


0.41+0.02 


0.055 


5.0 


G-14 


2.18±0.13 


3100 


7.30+0.07 


0.0169+0.0010 


0.452+0.014 


0.09 


6.4 


G-102 


1.07±0.17 


2900 


6.33+0.08 


0.031+0.002 


0.70+0.03 


0.07 


2.0 


CrAPMS 3B 


0.33±0.15 


2700 


6.14+0.07 


0.028+0.002 


0.77+0.03 


0.04 


1.1 


IRAC-CrA 3 


1.73±0.25 


2700 


8.48+0.13 


0.0033+0.0004 


0.261+0.015 


0.03 


9.0 



References. Lopez Marti et al. J2005I and |in preparation [l; Sicilia-Aguilar et al. ( i2008t . Neuhauser & Forbrich ((2008); Meyer &. Wilking i2009b 
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Fig. 4. Histograms of S values at ditferent temperature ranges in the Corona Australis star-forming region. The dashed lines indicate the expected 
contamination by field main sequence stars, and the hardly visible dotted lines at the bottom of each panel, the contamination by background cool 
giants. The presence of members of the star forming region is evident in those ranges of 5 where contamination is expected to be negligible. The 
hashed histograms show the contribution of objects from previous studies to these excesses. 



3 Myr and a distance of 130 pc, according to the Baraffe et al. 
(flW) tracks)^ 

Furthermore, the low density of objects in the young locus 
together with the clear separation from the rest of sources in 



' As a matter of fact, a couple of known, highly extincted cluster 
members in this temperature range are actually present in our survey 
(namely CrA 466 and G-87, both with Ay > 8 mag). However, our fits 
to their SEDs are not reliable, as they yield too high effective tempera- 
tures (around 4400 K), probably due to binarity or variability. 



our survey suggest that contamination from field stars must be 
very low in this area of the diagram. Most of the remaining ob- 
jects are placed well below the 20 Myr isochrone, indicating that 
they are foreground dwarfs characterized by high values of S . 
Background contamination would be caused by cool giants in 
the same range of effective temperatures as our objects of inter- 
est, but this contamination is expected to be very low due to the 
high galactic latitude of Corona Australis (b ~ 18°), and to the 
dark cloud itself. 
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Table 2. Optical and near-infrared photometry of our new candidate members 



Name 

CrAJ1901 11.6-364532.0 
CrAJ190139. 1-370016.8 
CrA J190151.7-371048.4 
CrAJ190207.6-371 156.5 
CrA J190239.2-365311.1 



R 

18.067±0.004 
21.07±0.06 
16.451±0.002 
19.615±0.013 
17.928±0.005 



15.779+0.005 
18.023+0.006 
14.500+0.003 
17.022+0.003 
15.654+0.005 



J 

13.375+0.026 
14.257+0.029 
12.326+0.030 
14.066+0.029 
13.190+0.024 



H 

12.496+0.028 
13.045+0.033 
11.423+0.033 
13.093+0.029 
12.471+0.022 



Ks 

12.129+0.026 
12.296+0.026 
11.047+0.030 
12.541+0.023 
12.060+0.023 



Table 3. SED fitting results for our new candidate members of Corona Australis 



Name 


Av (mag) 


Teff (K) 


5 


L(Lo) 


R{Re) 


M(Mo) 


age (Myr) 


CrA J1901 11.6-364532.0 


1.78+0.13 


3000 


7.20+0.06 


0.0161+0.0009 


0.471+0.013 


0.07 


4.5 


CrA J190139.1-370016.8 


6.20+0.14 


3200 


7.13+0.07 


0.0224+0.0014 


0.487+0.016 


0.13 


8.0 


CrA J190151.7-3710484 


1.96+0.09 


3400 


6.58+0.05 


0.047+0.002 


0.627+0.013 


0.3 


14.3 


CrA J190207.6-371 156.5 


3.62+0.13 


3100 


7.54+0.06 


0.0135+0.0007 


0.404+0.011 


0.1 


10.0 


CrA J190239.2-365311.1 


1.3+0.2 


2900 


7.08+0.09 


0.0157+0.0013 


0.50+0.02 


0.05 


3.2 



Following ICSL09 1 the foreground contamination was esti- 
mated by computing the amount of main-sequence stars ex- 
pected in our surveyed area for different temperature bins ac- 
cording to the local initial mass function (IMF) by Chabrier 



(1200 3'). For T,,// < 3400 K, the number of these objects ex- 
pected to be found in the young locus in Fig.[3]is three (3). The 
level of background contamination was estimated with the tem- 
perature and radii for cool giants provided by Fluks (1199 8^ and 
the volume density law from Wainscoat et al. ( 119921 ). From this 
estimation only one (1) background contaminant is expected for 
Teff < 3400 K. Because we count 23 objects in the area of in- 
terest of the (Teff, S) diagram, the amount of expected contam- 
ination is about 17%. 



3.3. Object selection 

The procedure to identify candidate cloud members is the same 
as that outlined in ICSL091 and is illustrated in Fig. [4] For four 
intervals of effective temperature we plotted the expected S- 
parameter distribution of foreground and background stars. As 
seen in Fig.|4] the number of foreground stars (dashed lines) de- 
creases with decreasing values of S , and with decreasing Teff. 
On the other hand, the number of background stars (dotted lines) 
is negligible in all but the highest temperatures. 

In general, our histograms are close to the predicted numbers 
of contaminating stars, although incompleteness due to the lim- 
iting magnitudes affects the bins with the highest values of S . 
In those ranges of S where the contamination is expected to be 
negligible, we interpret an excess of sources as objects that are 
probable members of the star-forming region. Based on the in- 
spection of Fig. |4]and to minimize the number of contaminants, 
the following criteria were used to select candidate members of 
the Corona Australis star-forming region: 

3100 < r,// < 3400 K, 5<8; 

2800 < r,y/ < 3100 K, S <9; 

2500 < Teff < 2800 K, S < 10; 

2200 < Teff < 2500 K, S < 10. 

In this way, 4, 13, 5 and 3 candidates respectively were se- 
lected in each temperature bin, amounting to 22 in total. The 
objects were then visually inspected to reject bad detections ei- 
ther due to bad pixels (two objects) or to source confusion in the 
cross-match (one object). Sources for which the fits derived ages 
older than expected for Corona Australis members were also re- 
jected as likely contaminants; concretely, we discarded two ob- 



jects whose derived ages were older than 100 Myr. For another 
source, the SED fitting suggests an age of 18 Myr, which we con- 
sider to be slightly too old to be a member of Corona Australis. 
The residual in the best fit for this latter source is relatively high 
iX^ ~ 0.4) due to its faintness (and hence large photometric er- 
rors) in the R-band (R > 20 mag). Because we do not have any 
additional information to support membership for this object, we 
chose to be conservative and reject it as a cluster member as well. 

This left us with 19 objects in the list of selected candidate 
members of Corona Australis. Thirteen (13) of them had been 
presented in previous works as candidate members of the region, 
and eight of these candidates had been confirmed as members 
through optical or near-infrared spectroscopy (see references in 
Table 3T|i. Another object (IRAC-CrA 3) had been selected by 
us as a candidate class 11 member of Corona Australis in a 
parallel study based on Spitzer photometry (Lopez Marti et al. 
in pr eparatio n! . Our fitting results for these fourteen objects are 
listed in Table 3.1 They have effective temperatures between 
2700 and 3400 K and masses betweeen 0.03 and 0.4Mq. Only 
two objects have masses above O.IMq. The optical extinction 
varies from to 2.2 mag. 

Together with the previously known candidates, our list still 
contains five objects that are likely members of the Corona 
Australis region, all of them with ages younger than 15 Myr ac- 
cording to the SED fitting, and none with reported detections 
in previous works. They have effective temperatures and masses 
within the same range as the previously known objects, but tend 
to be somewhat more extincted, with Ay values in the range 1- 
4 mag and one object with Av > 6 mag. The optical and near- 
infrared photometry for these new candidate members is given 
in Table|3.2| and their parameters from the SED fitting are listed 



3.2 



in Table 

Note that (excluding the two objects rejected due to pho- 
tometry problems and the one cross-matching misidentification) 
we classified three sources as likely contaminants, which agrees 
very well with our estimation from Sect. |3.2| On the other hand 
it must be remarked that we are missing members of the region 
whose S -parameter values are in the ranges affected by contami- 
nation in the histograms of Fig.|4] Indeed, some excess above the 
expected foreground contamination is observed in some bins to 
the left of our S cutoff value for Teff <3100 K (upper panels of 
Fig.|4]i, where we would expect to find more cluster members. In 
addition, as explained in Sect. |3.1| wrong fits would be obtained 
for objects affected by variability, binarity, strong accretion or 
infrared excess, and these sources would be wrongly placed in 
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Fig. 5. Same as Fig. [b] for the previous candidate members (trian- 
gles) and new candidate members from this work (circles) selected in 
Sect. [331 



the (Teff,S) diagram of Fig. [3] As a consequence, cloud mem- 
bers with primordial disks may be slightly underrepresented in 
our selected sample. 

In particular, other candidate members from the studies of 
Fernandez & Comeron (120011 1 object), Lopez Marti et al. 
( 120051 3 objec ts), Fo rbrich et al. ( 12007 1 2 objects) and Sicilia- 
Aguilar et al. ( 120081 2 objects) are detected in our survey with 
photometry in all RIJHKs bands, but are not selected with our 
method. Hence, they are not listed in Table 3.1 nor in Fig. [T| (but 
they are included in Figs. [3] and |4]l. For the previously known 
candidate members present in our survey, our success rate in re- 
covering them is thus about 65%, or 72% considering only ob- 
jects whose membership had been confirmed by spectroscopy; 
these percentages agree well with the 30% of members with poor 
fits found in the Lupus clouds ( CSL091 ). For the same reasons 
there could be some yet unknown members of CrA which are 
not selected by our criteria. 



4. Discussion 

4.1. Membership of the candidates 

We now discuss the membership status of our five objects. 
Figure |5] is a {Teff, S) diagram analogous to Fig.[3j but show- 
ing only the location of the members and candidate members of 
Corona Australis selected in this work. 

The spatial distribution of our sources is shown in Fig . [T| 
Both the previously known candidate members from Table 3.1 
and the new candidate members from Table l3.2l are located in the 
northern WFI field; the only exception is the candidate member 
CrA 205 from Lopez Marti et al. (j2005t . which is placed in the 
southern field. This distribution supports the membership of our 
new candidates to the Coronet cluster around the intermediate- 
mass star R CrA. 

Two of our candidates, CrA J190111.6-364532.0 and 
CrA J190151.7-371048.4, have a USNO-B counterpart within 
2". We compare their proper motions as reported in this catalog 
with those of previously known members of the Coronet clus- 
ter (see Table |4.l| . We considered only USNO-B measurements 
with a high probability estimator (>0.7) and measurements in at 
least three epochs (Monet et al. 120031) . The proper motion data 
are summarized in Table 14. l] 
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Fig. 6. Upper panel: p.^ versus cos^ from the USNO-B 1 catalogue 
for Coronet candidate members from previous works (triangles) and our 
new candidate members (circles). Lower panel: Spatial location of the 
same sources. The lines indicate their expected displacement in 10^ yr. 



With only 14 objects in total it is not possible to perform 
a statistically meaningful membership analysis based on proper 
motion. The upper panel of Fig. |6] shows versus cos 6 for 
all sources. In the lower panel we have plotted the current spa- 
tial locations of the sources and their expected displacement in 
10^ yr. While we cannot draw any firm conclusion on their mem- 
bership, this analysis shows that at least one of the candidates, 
CrA J190111. 6-364532.0, has a proper motion consistent with 
other members. For CrA J190151. 7-371048.4 there is at least 
some marginal agreement (within the errors) in its proper mo- 
tion components with respect to those of the majority of clus- 
ter members plotted in Fig. |6] Inspection of the optical images 
shows that this object has two close visual companions, located 
about 3" to the NE and 2" to the NW respectively, a fact that 
may be affecting the proper motion measurements. 

It is remarkable that one source, namely CrA 4109, displays 
clearly larger proper motion than the rest of probable cluster 
members plotted in Fig. |6] with /i ^ 109 mas/yr. The direc- 
tion of its expected displacement is also significantly diff'erent. 
Given the low object density in our surveyed region, source con- 
fusion seems unlikely. Moreover, the USNO-B R and / photom- 
etry (16.07 and 13.98 mag, respectively) is fully consistent with 
the WFI photometry for this object (R = 15.99 and / = 13.98; 
Lopez Marti et al. 2005). Hence we are quite confident that our 
identification of CrA 4109 with this proper motion source is cor- 
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Table 4. Proper motion information for CrA members and candidate members 



Name 


USNO-Bl 


yu„ cos S (mas/yr) 


fis (mas/yr) 


P 


Nobs 


CrA 205 


0526-0958780 


38±6 


-52±16 


9 


4 


CrA 466 


0530-0881950 


24+54 


-64±2 


7 


3 


CrA 468 


0529-0905310 


20±9 


-34±9 


9 


4 


CrA 4107 


0532-0801755 


8±1 


-22±6 


9 


4 


CrA 4108 


0532-0800872 


6±10 


-24±4 


9 


4 


CrA 4109 


0532-0800993 


108±5 


14±5 


9 


4 


CrA 41 10 


0530-0881945 


28±2 


-48±5 


9 


4 


CrAPMS 3AB 


0530-0882190 


54±20 


-28±7 


9 


4 


[FP2007] J190149.35-370028.6 


0529-0905310 


20±9 


-34±9 


9 


4 


G-1 


0530-0882197 


22±7 


-24±16 


9 


4 


G-14 


0529-0905404 


-8+1 


-24±22 


7 


4 


G-1 02 


0529-0905023 


12±8 


-44±12 


9 


4 


RXJ1901. 1-3648 


0532-0800238 


4±23 


-40±14 


9 


4 


CrA J 190 15 1.7-371048.4 


0528-0928314 


30±5 


-6±12 


9 


4 


CrA J1901 11.6-364532.0 


0532-0800263 


10±12 


-26±8 


9 


4 



Refere nces. Neuhauser et al. (12000 ): Lopez Marti et al. J2005L Forbrich & Preibisch ( I2007t : Sicilia-Aguilar et al. ( I2008t : Neuhauser & Forbrich 
(gDGD; this work 



rect. On the other hand, although there is no definitive confir- 
mation of its membership to Corona Australis, there are several 
indications in the literature (Ha emission, mid-infrared photom- 
etry and spectroscopy) suggestive of CrA 4109 being a young 
object (Lopez Marti et al. 2005; SiciHa-Aguilar et al. 120081) . 

In Fig. [T] CrA 4109 is located on the outskirts of the cloud 
core in an area of relatively low extinction. A possibility would 
thus be that this object was ejected from its birth site due to dy- 
namical interactions with other cluster members, maybe within 
a multiple system. We note that two more sources (CrA 4107 
and CrA 4108) are located in the same area as CrA 4109, but 
their proper motions agree very well with those of the majority 
of cluster members. 



4.2. Disks? 

Given that our candidates from Table 13. 2] had not been detected 
in previous surveys (the vast majority of them looking for sig- 
natures of accretion or disks), they are expected to be disk- 
less, or to have at most optically thin disks. To check this, we 
cross-matched our optical-NIR catalog with the Spitzer IRAC 
and MIPS 24 /vm observations from the study reported in Lopez 
Marti et al. fin pre paration ; see that work for details on the re- 
duction and analysis of these data). 

The mid-infrared photometry for our new candidates is sum- 
marized in Table |4.2| Four of them have [3.6] - [4.5] and 
[5.8] - [8.0] colors close to zero, as expected for class III sources 
(e.g. Allen et al. l2004l Hartm ann et al. l2005t . With the cri- 
terion proposed by Lada et al. ('20061) to make a classification 
of the disks based on the value of the SED slope ajuAc in the 
3.6-8.0 yum range, we find that two sources, CrA J190239.2- 
365311.1 and CrA J190207.6-371 156.5, are classified as disk- 
less (a,RAc < -2.56), while CrA J190151.7-371048.4 and 
CrA J190139. 1-370016.8 could be surrounded by anemic (or 
thin) disks (-2.56 < uirac < -1-80). We remark though that 
the slope values of the last two objects (-2.556 and -2.397, re- 
spectively) are very close to the defined limiting value between 
anemic disks and no disks. Furthermore these limits are not ab- 
solute, because they depend slightly on spectral type and (more 
strongly) on photometric errors (Lada et al. 2006 ; Hernandez 
et al. 120081 ). Given the dispersion on the slope of the field 
(thus, diskless) objects in our survey (cr ~ 0.1), we consider 



CrA J190151.7-371048.4 and CrA J190139.1-370016.8 as disk- 
less. The remaining source (CrA J1901 1 1.6-364532.0) could not 
be classified according to these schemes, because it was not ob- 
served at 4.5 and 8.0 fim. 

The SEDs of our candidates are shown in Fig.|7]compared to 
the best-fitting model photosphere for each object. We note that 
the error bars are smaller than the symbol size in the plots due 
to the low flux uncertainties. From visual inspection, the four 
sources classified as diskless display nearly photospheric SEDs 
in the whole range in which they are detected (up to 8 jum). 

The source CrA J190111.6-364532.0 is the only one with a 
counterpart at 24 jum. As seen in Fig.jTj its SED has the charac- 
teristic shape of a star that has cleared its inner disk, showing a 
clear excess over the photosphere only after 10 jum — a so-called 
"transitional disk" (e.g. Calvet & D'Alessio '2001', Calvet et al. 
[2005). This source would have a mass of about O.O7M0 accord- 
ing to the SED fitting, which would make it one of the lowest- 
mass objects repotted to possess a transitional disk to date. 

4.3. Disk fraction and transition disks 



Among the previously known objects selected in Sect. 3.3 and 
listed in Table |3.l| there are five class II sources (G-14, CrA 465, 
CrA 4107, CrA 4110 and IRAC-CrA 3), four class III sources 
(CrA 453, CrA 468, CrA 4108 and CrA 411 if] and three tran- 
sitional objects (G-1 02, Cr A 205 and CrA 4109) according to 
Sicilia-Aguilar et al. (120081 ) and our own analysis (Lopez Marti 
et al. |in preparation i. Note that by "transitional object" we mean 
here a star or brown dwarf harboring a disk with an inner hole, as 
suggested by a 24 fim excess larger than at 8 //m or shorter wave- 
lengths. No classification is available for CrAPMS 3B, because 
this is the low-mass companion of a PMS star and is unresolved 
in the IRAC images, nor for CrA 452, which lies outside the area 
covered by the IRAC observations. 



- CrA 4111 was classified as a transitional object (in the sense in 
which we are using the term here) by Sicilia-Aguilar et al. (2008 ), who 
report the detection of a 24 yum excess for this source. Our own analysis 
of the Spitzer data, however, does not confirm this result: We do not find 
any counterpart for CrA 4111 at 24 yum. Visual inspection of the MIPS 
24 yum mosaic shows that this object is located in an area of nebulosity, 
which could have yielded a spurious detection. 
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Fig. 7. Spectral energy distributions of our new candidate members of Corona Australis (red dots), from the optical to the mid-infrared, compared 
to their corresponding best-fitting photospheres according to the SED analysis (see model parameters in Table|3.2|(. 



If we now add our five new candidates to the sample, the final 
census amounts to five sources with optically thick disks, eight 
sources without disks, and four sources with transitional disks 
for a total of 17 classified objects selected following the same 
criteria. Taking both primordial and transitional disks together, 
the disk fraction in our sample amounts to 53±18%, lower than 



the 70% reported by Sicilia-Aguilar et al. ( I2008I I, due to the in- 
crease of diskless objects. 

The primordial disk fraction in our sample is 30± 13%, lower 
than the estimations of Haisch et al. (2001 ) and Mamajek (20()9]l 
for a 3 Myr cluster (the usually quoted age for the Coronet clus- 
ter), but not far from their given value of about 50+10% when 
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Table 5. Spitzer photometry of our new candidate members 



Name 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24.0] 


aiRAc 


Class 


Disk" 


CrAJ1901 11.6-364532.0 


11.653+0.004 




11.375+0.015 




8.23+0.075 








CrAJ190139.1-370016.8 


11.653±0.005 


11.48±0.006 


11.415+0.016 


11.392+0.026 


n.d. 


-2.556 


III 


A/N 


CrA J190151.7-371048.4 


10.553±0.003 


10.44±0.003 


10.265+0.007 


10.182+0.008 


n.d. 


-2.397 


III 


A/N 


CrAJ190207.6-371 156.5 


11.893±0.004 


11.86±0.006 


11.795+0.017 


11.902+0.03 


n.d. 


-2.836 


III 


N 


CrA J190239.2-365311.1 


11.673±0.004 


11.54±0.005 


11.455+0.015 


11.492+0.018 


n.d. 


-2.633 


III 


N 



Notes. n.d.=not detected; blank=not observed 
A=thin disk; N=no disk (see discussion in text) 



the errors are considered. We note, though, that a relatively large 
age spread is observed among the objects in our sample, ac- 
cording to the SED fitting results. The derived primordial disk 
fraction is only slightly higher, within the errors, than that ob- 
tained by Haisch et al. (2001) for the 5 Myr NGC 2362 cluster 
(12 ± 4%), which agrees quite well with the median age of our 
sample (5 IVIyr). However, this result should not be interpreted 
as a new estimation of the cluster age, for several reasons. The 
first one is of course the low number of objects in our sample and 
its limited mass range. Second, we caution that individual ages 
derived from theoretical SED fitting are affected by the many 
uncertainties still present in the models and may even depend on 
the accretion history of the objects (e.g. Baraffe et al. l2002|[2009l 
Mayne & Naylor 2008 ). In addition, as explained in Sect. [3] the 
objects with strongest accretion (thus probably the youngest) are 
likely to be missed by our selection, which produces a slightly 
older sample than the mean age of the cluster. The standard devi- 
ation of our sample age (3.9) accounts well for these uncertain- 
ties. 



On the other hand, the ratio of transitional to primordial disks 
(45±22%) agrees quite well with the 50% value reported by 
Sicilia-Aguilar et al. (but see also Ercolano et al. 12009b . As al- 
ready noted by these authors, this ratio is higher than observed in 
other regions of similar age (around 5-10%; see e.g. IVTerin et al., 
|in preparation^ for a discussion of the transitional disk ratios in 
the c2d clouds). However, the majority of these studies targeted 
higher mass stars than considered in Sicilia-Aguilar et al. ( 12008b 
and the present study. The larger number of transitional disks 
among low-mass stars and brown dwarfs with respect to solar- 
type stars could then be interpreted as a hint of a faster disk evo- 
lution in the former, as suggested by previous works (e.g. Sterzik 
et al. 120041 Bouy et al. 120071) . It also suggests that the "transi- 
tional" stage is not necessarily a short-lived step in the evolution 
towards a protoplanetary or debris disk, at least for very low- 
mass objects. However, our sample is too small to derive any 
statistically meaningful conclusion on this issue. 



We remark that the inclusion of the two unclassified objects 
(CrAPMS 3B and CrA 452) in this census will not alter the mea- 
sured disk fraction by more than 5% in either sense (higher or 
lower disk fraction), nor will it change the ratio of transitional to 
primordial disks by more than 10%, and hence will not substan- 
tially affect these conclusions. On the other hand, as explained in 
Sect. 



3.1 



cloud members with primordial disks may be slightly 
underrepresented in our selected sample. However, even in the 
most pessimistic case (that would be, if the whole 30% of badly 
fitted sources estimated in Sect. 3.1 were surrounded with pri- 
mordial disks), the disk fractions would stiU be within the un- 
certainties of our estimations. 



5. Conclusions 

We performed an analysis of optical and near-infrared data of the 
Corona Australis star-forming region based on the S -parameter 
formalism. In our surveyed area of ~ 0.64 deg^, we identified 
fourteen previous candidate members and five new candidate 
members of this dark cloud. 

The new candidates have estimated effective temperatures 
between 2900 and 3400 K, corresponding to masses between 
0.05 and 0.13Mq. They are thus very low-mass stars and massive 
brown dwarfs. Their ages span between 3 and 15 Myr, which is 
consistent with the reported age spread in Corona Australis (~1- 
10 Myr), given the uncertainties in the models, in the distance 
to the cloud and in the SED fitting procedure. The membership 
of these objects to the star forming region is further supported 
by their spatial distribution and, when available, proper motion 
information. The source for which membership is more uncer- 
tain is CrA J190151.7-371048.4, the oldest star in our sample 
according to the SED fitting results. 

The SEDs of four of our new candidates are nearly photo- 
spheric. The exception is CrA J1901 1 1.6-364532.0: This source 
displays excess at 24 /im, which is suggestive of a transition disk 
with an inner hole. With an estimated mass of 0.07 Mq, this is one 
of the lowest-mass objects reported to possess such a disk. 

We calculated the disk fraction of the Corona Australis pop- 
ulation selected with our method, which is 50%. This value 
is lower than the one reported in a previous study by Sicilia- 
Aguilar et al. (12008b . The ratio of transitional to primordial disks 
(45%) though agrees well with the fraction of 50% reported by 
Sicilia-Aguilar et al. (2008 ) and is remarkably higher than the 
value measured in other clusters of similar age. This suggests 
that transitional disks around brown dwarfs may have longer 
lifetimes than around low-mass stars. However, this impression 
should be confirmed with a larger, more statistically meaningful 
sample. 

The results from this work stress the need to properly char- 
acterize the diskless population of a region to derive meaningful 
disk fractions. This is especially important to understand the de- 
pendence of the disk fractions with the mass of the central ob- 
jects, in particular for the very low-mass population. Even in a 
relatively small survey like the one reported here and in a well- 
studied and relatively low-density region like Corona Australis, 
25% of the total number of cloud members selected with our 
method had not been detected in previous studies based on accre- 
tion, disk or activity signatures. An analogous study of an even- 
tual larger-scale survey covering most of the dark cloud would 
probably identify a significant number of young very low-mass 
objects belonging to this star-forming region. 
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